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I. INTRODUCTION 


In the Standard Model (SM), the weak charged-current 
transitions mix quarks of different generations, which 
is encoded in the unitary Cabibbo-Kobayaslii-Maskawa 
(CKM) matrix [lj 2|. In the case of three generations of 
quarks, the physical content of this matrix reduces to four 
real parameters, among which one phase, the only source 
of CP violation in the SM (neglecting CP-violating ef¬ 
fects induced by the strong-interaction 0-term or neutrino 
masses): 


A 2 
p + ifj 


\V US \ 2 

\v ud \ 2 + \v us \ 2 ’ 

V udV* b 
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Wcb \ 2 
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( 1 ) 


One can exploit the unitarity of the CKM matrix to de¬ 
termine all its elements (and when needed, to obtain their 
Wolfenstein expansion in powers of A) [3HS]. 

Extracting information on these parameters from data 
is a challenge for both experimentalists and theorists, 
since the SM depends on a large set of parameters 
which are not predicted within its framework, and must 


be determined experimentally. An additional difficulty 
stems from the presence of the strong interaction bind¬ 
ing quarks into hadrons, which is responsible for most 
of the theoretical uncertainties discussed when determin¬ 
ing the CKM matrix parameters. The CKMfitter group 
aims at this goal by combining a large set of constraints 
from flavour physics, using a standard y 2 -like frequen¬ 
tist approach, in addition to a specific (Rfit) scheme to 
treat theoretical uncertainties 0E| (see refs. UHTO] for 
alternative approaches in this context). 

As will be illustrated below, the SM global fit has 
reached a remarkable accuracy from both the experimen¬ 
tal and theoretical points of view. In this context, and 
following a long history of flavour as a probe for “New 
Physics” (existence of the charm quark, bounds on the 
top quark mass...), one can also use flavour observables 
to constrain models of New Physics (NP), either in a par¬ 
ticular scenario or with a rather generic scope. We will 
follow the second avenue, providing results for generic 
New Physics in A F = 2 and updating ref. m- 

The results presented here correspond to the most re¬ 
cent update performed by the CKMfitter collaboration, 
including results obtained until the CKM 2014 workshop 
in Vienna [5|. 
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TABLE I. Constraints used for the global fit, and the main inputs involved (more information can be found in ref. [6]). When 
two errors are quoted, the first one is statistical, the second one systematic. The lattice inputs are our own averages obtained 
as described in the text. 


CKM 

Process 

Observables 

Theoretical inputs 

\V u d\ 

0 + —> 0 + transitions 

| Vud | nucl 

= 

0.97425 ± 0 ± 0.00022 

m 

Nuclear matrix elements 

\V US \ 

K —> 7 rtv 

Ws SL/+ _> 7 r ( 0 ) 


0.21664 ±0.00048 

m 

if^(o) 

= 0.9641 ±0.0015 ±0.045 


K —> ev e 

B[K -> ev e ) 


(1.581 ±0.008)- 10“ s 

El 

/if 

= 155.2 ± 0.2 ± 0.6 MeV 


K ->• pv^ 

B(K -»■ pvJ 

= 

0.6355 ±0.0011 

m 




T — > I\V T 

B(t ->• Kv T ) 

= 

(0.6955 ± 0.0096) ■ 10 " 2 

m 



IW.I 

Wul\ 

K —> pv/ir —> pv 

B(K -> pv 

B( TV —> pV M ) 

= 

1.3365 ±0.0032 

na 

SkIU 

= 1.1942 ±0.0009 ±0.0030 


t — r Kv/t —¥ nv 

B(t Kv t ) 
B(t —¥ nv T ) 

= 

(6.43 ±0.09) ■ 10 ” 2 

m 



\Vcd\ 

vN 

| Vcd | uN 

= 

0.230 ±0.011 

na 




D —» pv 

B(D -)• pv) 

= 

(3.74 ±0.17) ■ 10 " 4 

m 

f Ds / f D 

= 1.201 ±0.004 ±0.010 


D —¥ irlv 

\V cd \f?^(0) 

= 

0.148 ±0.004 

m 

0) 

= 0.666 ± 0.020 ± 0.048 

\Vc\ 

W —¥ cs 

| Vcs | W —>cs 


0.94/°; 22 ±0.13 

na 




D s —> TV 

B(D S -»■ tv) 

= 

(5.55 ±0.24) ■ 10 ” 2 

m 

fDs 

= 245.3 ± 0.5 ± 4.5 MeV 


D s —> pv 

B(D S -¥ pv^) 

= 

(5.57 ±0.24) ■ 10 ” 3 

na 




D ->• Kiv 

\v„\f?^ K (p) 

= 

0.712 ±0.007 

mm 

f?^ K ( 0) 

= 0.747 ±0.011 ±0.034 

\v ub \ 

semileptonic decays 

IWblsL 

= 

(3.70 ±0.12 ±0.26) • 10 -3 

na 

form factors, shape functions 


B —¥ TV 

B(B tv) 

= 

(1.08 ± 0 . 21 ) ■ 10' 4 

nana 

/s a //s 

= 1.205 ±0.004 ±0.007 

\Vcb\ 

semileptonic decays 

| Vcb ISL 

= 

(41.00 ±0.33 ±0.74) • 10 “ 3 

na 

form factors, OPE matrix elements 

a 

B —¥ 7T7T, p%, PP 

branching ratios, CP asymmetries 

na 

isospin symmetry 

p 

B -¥ (cc)K 

sin(2/3) [cS] 

= 

0.682 ±0.019 

na 


7 

B -¥ 

inputs for the 3 methods 

na 

GGSZ, GLW, ADS methods 


B a ->■ J/lj>(KK , 7T7T) 


= 

-0.015 ±0.035 

na 


v: q v w 

A m d 

A m d 

= 

0.510 ±0.003 ps " 1 

na 

B B s/B Bd 

= 1.023 ±0.013 ±0.014 


A m a 

Arris 

= 

17.757 ±0.021 ps " 1 

na 

b Bs 

= 1.320 ±0.017 ±0.030 


B s —y nn 

B(B a —> pp) 

= 

( 2 . 8 ±g; 5 ) ■ 10" 9 

na 

fa a 

= 225.6 ±1.1 ±5.4 MeV 

V t * d v ts 

£K 

M 

= 

(2.228 ± 0 . 011 ) • 10“ 3 

na 

Bk 

= 0.7615 ±0.0027 ±0.0137 

v: d v C3 






Ke 

= 0.940 ±0.013 ±0.023 


II. INPUTS FOR THE SM GLOBAL FIT 

A. General discussion 

Not all the observables in flavour physics can be used 
as inputs to constrain the CKM matrix, due to limita¬ 
tions on our experimental and/or theoretical knowledge 
on these quantities. The list of inputs to the global fit 
is indicated in Table [j] they fulfill the double require¬ 
ment of a satisfying control of the attached theoretical 
uncertainties and a good experimental accuracy of their 
measurements. In addition, we only take as inputs the 
quantities that provide constraints on the CKM parame¬ 
ters A, A, p, fj. We will see below that not all parameters 
are equally relevant for the global fit. 

A major source of uncertainties in flavour analyses 
arises from matrix elements that encode the effects of the 
strong interaction in the nonperturbative regime, corre¬ 
sponding here to decay constants, form factors and bag 
parameters. We rely mainly on lattice QCD simulations 
for the determination of these quantities, as they provide 


well-established methods to compute these observables 
with a controlled accuracy. Some of the uncertainties 
have a clear statistical interpretation. Lattice simulations 
evaluate Green functions in a Euclidean metric expressed 
as path integrals using Monte Carlo methods, and their 
accuracy depends on the size of the sample of gauge con¬ 
figurations used for the computation. The remaining 
uncertainties are systematic: they are now dominant in 
most cases and they depend on the computational strate¬ 
gies chosen by competing lattice collaborations: discreti¬ 
sation methods used to describe gauge fields and fermions 
on a lattice, interpolating fields, parameters of the sim¬ 
ulations, such as the size of the (finite) volumes and lat¬ 
tice spacings, the masses of the quarks that can be sim¬ 
ulated, and the number of dynamical flavours included 
as sea quarks. These simulations must often be extrapo¬ 
lated to obtain physical quantities, relying in particular 
on effective theories such as chiral perturbation theory 
and heavy-quark effective theory which induce further 
systematics. 

The combination of lattice values is a critical point of 
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most global analyses of the flavour physics data, even 
though there is no universal definition of theoretical un¬ 
certainties (and hence how to combine them). Several ap¬ 
proaches have been proposed to perform such a combina¬ 
tion Gang, and we have also proposed our own scheme, 
systematic, reproducible and to some extent conserva¬ 
tive m ■ We have collected the relevant lattice calcula¬ 
tions of the meson decay constants, as well as the Bd , 
B s and I\ bag parameters, and the I\ —* ir, D — > tt and 
D —> K vector form factors at zero momentum trans¬ 
fer. We base our set of calculations on the latest FLAG 
(Flavour Lattice Averaging Group) report |I5], with the 
addition of new results published since that report was 
written [6]. We perform our averages considering values 
from lattice simulations with different numbers of dy¬ 
namical flavours (Nf = 2, 2 + 1, 2 + 1 -f 1). Even though 
the different collaborations attempt at assessing the cor¬ 
responding systematics in a careful way, one cannot ex¬ 
clude that such combinations are affected by further sys¬ 
tematics which unfortunately cannot be assessed easily. 
These lattice averages are the input parameters used in 
the fits presented in this paper. In the specific case of de¬ 
cay constants, the S'{7(3)-flavour breaking ratios /k/ f-n, 
IdJ/d: f b s IfBa are better determined than the individ¬ 
ual decay constants. We will therefore take these ratios 
as well as the strange-meson decay constants as reference 
quantities for our inputs. In the same spirit, it is more 
relevant to consider the predictions of the ratio K^l^n 
of the kaon and pion leptonic partial widths, as well as 
B(t — > Kv t )/B(t — » 7 tv t ) instead of individual branch¬ 
ing ratios. 

There are also other sources of theoretical uncertain¬ 
ties. This is the case for the inclusive and exclusive deter¬ 
minations of \V u b\ and \V c b\, which involve nonperturba- 
tive inputs of different natures. We use the latest HFAG 
results m for each of these determinations and com¬ 
bine inclusive and exclusive determinations following the 
same scheme as for the combination of lattice quantities. 
We also need theoretical inputs for heavy up-type quark 
masses, namely fh c (fh c ) = (1.286 ± 0.013 ± 0.040) GeV 
and fhtifht) = (165.95±0.35±0.64) GeV, as well as to the 
strong coupling constant a s (Mz ) = 0.1185 ± 0 ± 0.0006. 
We refer the reader to refs. PDH25] for a more detailed dis¬ 
cussion of each constraint, whereas the related hadronic 
inputs can be found in ref. [5]. 


B. Specific inputs 

A few specific inputs have changed recently and de¬ 
serve comment. 

Constraints on the CKM angle a are derived from 
the isospin analysis of the charmless B ±, ° —> (7r7r) ±,cl , 
B ± '° —>• (pp) ±fi and B° —► (p7r)° decays. Assuming the 
isospin symmetry and neglecting the electroweak penguin 
contributions, the amplitudes of the SU(2)-conjugated 
modes are constrained by triangular (or pentagonal) re¬ 
lations. The measured branching fractions and asymme- 



FIG. 1. Inputs for a used in the SM global fit. We show the 
world averages for mr, pn and pp separately. 

tries in the B ± -° —► (tttt ) 1 * 1 ’ 0 and B ± ’° —► (pp ) 1 * 1 ’ 0 modes 
and the bilinear form factors in the Dalitz analysis of 
the B° —► (p 7 r)° decays provide enough observables to 
simultaneously determine the weak phase /3 + 7 = 7 r — a 
together with the tree and penguin contributions to each 
mode. 

The world average constraint on a is so far dominated 
by the B ± ’° —>• (pp) ±,c) data, thanks to the low level of 
the penguin contribution to these modes, conducting to 
the 68.3% Confidence Level (CL) intervals : 

a(B -7 pp) = (89.9 tti)° U (0.1+g;|)°. (2) 

The recent update of the measured branching fraction of 
the B —» 7 T 0 7 r° decay, driven by the Belle experiment [M] , 
significantly improves the determination of a through 
the isospin analysis of the B ± -° —> ( 7 T 7 r) ± 0 modes. The 
68.3% CL intervals 

a(B 7T7r) = (95.0i?‘|)° U (175.0j£|)° U (135.5 ± 15)° 

(3) 

are obtained. Combining the experimental data for the 
7 T 7 T, pp and p 7 r decay modes, the world average 68.3% CL 
intervals 

%a = (87.7±i;i)° U (179.0±|;I)° (4) 

are obtained (Fig. [I]). The recent Belle update on 
B(B° —> improves the la a resolution by 0.5° with 

respect to the previous determination. 

For the constraint on 7 , we have considered recent re¬ 
sults from l?-meson decays to open-charm final states, 
B —> £)(*)/<"(*). In those decays, the interference between 
b —» cus and b —> ucs tree amplitudes gives access to 
the weak phase 7 . Several methods have been proposed, 
which can be grouped according to the choice of the final 
state. Recent results include the updated LHCb results 
for the charged B —>• DK decay, where D —)■ Kstt + tt ~, 
KsK + K~, using a 3 fb _1 data sample [25j and for the 
first time, several observables, including CP asymme¬ 
tries, for the B° —► DK*° decays, where D decays in 
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FIG. 2. Inputs for 7 used in the SM global fit (top). We show 
the world averages for the different methods, in the ( 7 ,S B ) 
(middle) and ( 7 , re) (bottom) planes. Shaded areas (dashed 
lines) enclose points with 1 — p < 68.3 % (95.45 %). 


FIG. 3. Constraint obtained for the r B parameter of the B —> 
DK*° mode. “LHCb ( DK * 0 )” (dot-dashed line) includes only 
data from DK*° whereas “LHCb” (dashed line) involves all 
channels (including DK *°). 


Parameter 

Value and uncertainties 

charged B 


r B (DK) 
5 b (DK ) 
r B (D*K) 
5 b (D*K) 
r B (DK*) 
5 b (DK *) 

0.0970t° 0 °° 0 6 63 

(125.4±™)° 

0 . 11918:819 

(-49118)° 

PI 1 Q’7+0-051 
u - ±0 ' _0.047 

(112112)° 

neutral B 


r B (DK*) 
5 b (DK *) 

U.ZiOU_q 052 
(336l^)° U (200lg°)° 


TABLE II. Confidence intervals for the main hadronic param¬ 
eters obtained from the combination of the relevant BaBar, 
Belle and LHCb observables measured in the charged and 
neutral B —> decays. 


tt~K + , K~I< + TT~n + [26]. Combining the experimental 
data for the decay modes, the world average 68.3% CL 
interval 


Iwa = (73.2 i?;g)° (5) 

is obtained (Fig. [2]), as well as the hadronic parameters 
(rs, the magnitude of the ratio of the amplitudes, and 
5b, the relative strong phase between the two ampli¬ 
tudes) summarized in Table |TT] Though the impact of 
the observables for the neutral B decay B - 7 - DK*° is 
small for the 7 measurement itself, it is worth noticing 
that the corresponding rs is now clearly measured away 
from 0, as rsi.DK* 0 ) = 0.236to;o52 (Fig.j3j). The recent 
measurement of LHCb with the B s —> u s K mode m 
has not been included in our 7 average. Though very 
promising while using only 1 fb^ 1 , we estimate its im¬ 
pact on the 7 error to be at the order of 0 . 1 °. 

Other quantities which have experienced recent im¬ 
provement are the branching ratios B(B q —>■ /r/i) with 
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FIG. 4. Constraints on the CKM coordinates with M = db, sb, ds,ct,ut,uc, from the global SM CKM-fit. Regions 

outside the coloured areas have 1 — p > 95.45 %. For the combined fit the yellow area inscribed by the contour line represents 
points with 1 — p < 95.45 %. The shaded area inside this region represents points with 1 — p < 68.3 %. 
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q = d, s. B(B S —> pp) have been observed and mea¬ 
sured both by CMS and LHCb (at 4.3er and 4.0cr respec¬ 
tively), leading to a rather accurate combination T5] . 
There are also interesting information already available 
on B(B d —► pp), even though the threshold for evidence 
has not been reached yet by either of the two experi¬ 
ments. On the theoretical side, new computations have 
been performed including NLO electroweak corrections 
and NNLO strong corrections settling down is¬ 

sues met by earlier calculations concerning the stability 
with respect to higher-order corrections. In our predic¬ 
tions, we include the residual uncertainty of 1.5% dis¬ 
cussed in ref. j2Sj. We will predict the value of the dilep- 
tonic branching ratios without time integration, which 
would induce a further increase of 0(AT s /r s ), more pre¬ 
cisely (1 + y s ) = 1.07 discussed in refs. pH1 - i55] . 


III. RESULTS OF THE SM GLOBAL FIT 


A. CKM parameters and Unitarity Triangles 


The current situation of the global fit in the (p, fj) plane 
is indicated in Fig. [4] Some comments are in order be¬ 
fore discussing the metrology of the parameters. There 
exists a unique preferred region defined by the entire set 
of observables under consideration in the global fit. This 
region is represented by the yellow surface inscribed by 
the red contour line for which the values of p and fj with 
a p-value such that 1 — p < 95.45 %. The goodness of the 
fit can be addressed in the simplified case where all the 
inputs uncertainties are taken as Gaussian, with a p -value 
found to be 66 % (i.e., 0.4 er; a more rigorous derivation 
of the p-value in the general case is beyond the scope of 
this article 05). One obtains the following values (at 
1 <t) for the 4 parameters describing the CKM matrix: 

A = 0.810 t°;°“ , A = 0.22548 , ( 6 ) 

P = 0-145 ±g; 8 S?, 7/ = 0.343™2- (7) 


The various constraints can be expressed in the unitarity 
triangles associated with the different mesons of interest, 
with angles defined independently of phase conventions: 


a dl d 2 = arg 


7 dl d 2 = arg 


Vt dl v ; d2 

Vink v : d2 

Vud,! Vud 2 
Vc dl v ; d2 


Pd ld2 = arg 


Vc dl v : d2 

v tdl v t * d2 


( 8 ) 


and similarly for the angles in the up sector: 


= arg 


7 uiu 2 = arg 


Vu lb v: 2b 

T/ T/* 

v u ld v U2 d 

Vu ld Vu 2d 


, p UlU2 = arg 


V^sV^s 

Vu lb v: 2b 


(9) 


One recovers the usual <j>i, <f> 2l <$>3 and a , /?, 7 (without 
subscripts) for the B d Unitarity Triangle (di = d, d 2 = b). 


B s —W 
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FIG. 5. Pulls for the SM global fit obtained by comparing the 
value of Xmin with and without including the measurement of 
the quantity. Notice that the different pulls are not necessarily 
independent. 


In the same general way the relative coordinates of the 
upper appex of each triangle are defined as 

- , ■- _ ^ dl V*d 2 

Pd\d 2 + ' l ' r ] dld 2 t r r r* ! 

V cd! V cd2 


- 1 •- ’“lu'Koa 

Pu 1 U 2 + ' lr lu 1 U 2 — ~ TP 7/7 ! ( 10 ) 

V UiS V U 2 S 

where again p+ifj = p d b + ifj d b refer to the B d system. In 
the B s case, (j) s can be defined as 2/3 s b- The corresponding 
triangles are shown in Fig. [4j in particular the ( sb ) where 
the constraint from <j) a is shown (but the corresponding 
label is not indicated). 


B. Comments and predictions 

As underlined above, the overall consistency seen 
among the constraints allows us to perform the metrol¬ 
ogy of the CKM parameters and to give predictions for 
any CKM-related observable within the SM. Let us add 
that the existence of a 1 —p < 95.45 % region in the (p, fj) 
plane is not equivalent to the statement that each individ¬ 
ual constraint lies in the global range of 1 — p < 95.45 %. 
Each comparison between the prediction issued from the 
fit and the corresponding measurement constitutes a null 
test of the SM hypothesis. 
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TABLE III. Comparison between prediction and measurement of some flavour observables in the SM. The first column describes 
the observables. The second and third columns give the measurement and the prediction from the global fit (not including 
the measurement of the quantity considered), respectively. The fourth column expresses the departure of the prediction to the 
measurement, when available. 

Observable Measurement Prediction Pull (a) 


Charged Leptonic Decays 


B(B+ 

->• t + v t ) 

(10.8 ±2.1) ■ 

10 -5 

mm 

(7.58 

+0.80 \ 
-0.59/ ' 

10 -5 

1.5 

B(B+ 

->• 

< 1.0- 10 

-6 

na 

(3.64 

+0.27 \ 
-0.38/ ' 

10- 7 

- 

B(D+ 

T + Vr) 

(5.55 ±0.24) 

10 -2 

m 

(5.19 

+0.02 \ 
-0.12) ' 

10 -2 

1.5 

B(D+ 


(5.57 ±0.24) 

10 -3 

m 

(5.31 

+0.02 \ 
-0.09/ ' 

10 -3 

1.1 

B{D+ 


(3.74 ±0.17) 

10~ 4 

na 

(3.91 

±0.11) 

■ 10“ 4 

0.6 


Neutral Leptonic B decays 


B(B? 

—» 

T + T ) 

- 



(6.92 

+0.41 \ 
—0.39/ 

10’ 

-7 

- 

B(B° S 

—> 

P + P~) 

(2-8 

■ 10" 9 

m 

(3.34 

+0.13\ 

—0.25/ 

10' 

-9 

1.0 

B(B° S 

->• 

e + e~) 

< 2.8 • 

10" 7 

na 

(7.64 

+0.46 \ 
-0.43J 

10' 

-14 

- 

B(B° d 

-* 

T + T~) 

< 4.1 • 

10" 3 

ca 

(2.05 

+0.13 \ 
-0.14/ 

■ 10 

-8 

- 

B(B° d 

—» 

P + P~) 

(3.6+1) 

■ 10” 10 

m 

(0.98 

+0.06 \ 
-0.07/ 

10' 

-10 

- 

B(B° d 

-1 

e + e~) 

< 8.3- 

10" 9 

na 

(2.29 

+0.14\ 

-0.16/ 

10' 

-15 

- 


B q — B q mixing observables 


AIL (ps- 1 ) 

0.081 ±0.006 

m 

0.120 

+0.043 

-0.045 

0.1 

a SL 

(1 ± 20) ■ 10~ 4 

na 

(-6.5 

±i:i ) • 10- 4 

0.3 

a SL 

(-48 ±48) • 10 -4 

na 

(0.29 

1 

O 

00 00 
0 0 
6 d 
+ 1 

1.0 

Asl 

(-47 ± 17) • 10" 4 

Ea 

(-3.4 

+ 1.0) . 10 -4 

2.7 

sin(2/3) 

0.682 ± 0.019 

m 

0.771 

+0.017 

-0.041 

1.7 

0s 

-0.015 ±0.035 

m 

-0.0365 

+0.0013 

-0.0012 

0.6 


Rare K decays 


B(K + —> n + i/v) 

(l-75±i;5g) ■ 10- 10 m 

(0.85 lo.ia) 

10- 10 

0.7 

B(K l — ► 71 °vv) 

- 

(0-28 tg' 94 ) 

10- 10 

- 


Some of the corresponding pulls are reported in Ta¬ 
ble |III| and shown in Fig. [5j showing that there is no 
sign of discrepancy with our set of inputs. One should 
also notice that some of the quantities included in our 
fit have only a limited impact on the outcome. This is 
for instance the case for quantities where the measure¬ 
ment is compatible, but less precise than the SM predic¬ 
tion, like 4> s , B(B S —> p/x), or semileptonic and leptonic 
-D( s ) decays. In Table Till we also include observables 


that were not used as input constraints, either because 
they are not measured at a sufficient accuracy yet, e.g., 
B(Bd —>• ( + £~), or because the control on the theoretical 
uncertainties remains under discussion, e.g., AF S . The 
corresponding predictions can then be directly compared 
with their experimental measurements (when they are 
available). 


Before moving to specific observables and correlations, 
we briefly discuss the lasting discrepancy between deter¬ 
minations of \V ub \ and \V cb \ using exclusive and inclu¬ 
sive semileptonic decays. As indicated previously, the 
global SM fit is based on an average for the two ma¬ 
trix elements, taking into account the differences between 
statistic and systematic errors. In Fig. [6j we illustrate 
the results obtained by considering only exclusive (top) 


or inclusive (bottom) determinations for both \V ub \ and 
\V cb \. As expected, the constraint from \e K \ changes sig¬ 
nificantly due to the variation in \V cb \, whereas the \V ub \ 
constraint from B —► tv is found in better agreement 
with the inclusive input than the exclusive one. An ad¬ 
ditional interesting feature in the inclusive case is the 
appearance of a partial ring from the combined contri¬ 
bution of A rrid and Am s . This feature appeared already 
in the SM fit for Summer 2012, and can be explained 
by the fact that this constraint combines constraints on 
\V td \ 2 = A 2 A 6 [(l-p) 2 ± + ±0(A 4 )] (yielding a ring in the 
(p, ff) plane) and on \V ts \ 2 = A 2 A 4 [1 —A 2 (l —2p) 2 ±0(A 4 )] 
(cutting too large values of p). The overall agreement be¬ 
tween the various constraints remains excellent in both 
inclusive and exclusive fits, with very little variation in 
the global p- value at the best-fit point and the confidence 
interval for the four Wolfestein parameters with respect 
to the global SM fit obtained from an average of inclusive 
and exclusive values for \V ub \ and |V^f>|. 

























p 


FIG. 6 . Constraints on the CKM (p, fj) coordinates from the 
global SM CKM-fit using only exclusive (top) and inclusive 
(bottom) determinations of \V u b\ and \V c b\ from semileptonic 
decays as inputs. 


C. Specific observables and correlations 

We focus now on some specific observables and their 
correlations. A first example is given by the two- 
dimensional comparison for B(B —» tv) and sin 2/3 in 
Fig.0 showing that the discrepancy that used to affect 
the SM global fit [2D] has now been alleviated to a large 
extent (remaining only at 1 . 6 cr). As discussed in ref. 120] . 
this discrepancy had an experimental origin, and it has 
been reduced thanks to the addition of new data (the re¬ 
maining discrepancy is driven by the larger BaBar result 



sin 2p 


FIG. 7. Prediction on B(B —» tv) and sin 2,0 coming from 
the global fit (without the corresponding inputs) compared 
to current experimental information (cross). Regions outside 
the coloured areas are excluded at 1 — p > 95.45 %. 

compared to Belle measurement). 

One can also consider B(Bd, s —> /j,(i) as shown in Fig. [ 8 } 
showing the confidence contours from the combination 
of CMS and LHCb [IS]. One notices that NLO and 
NNLO predictions follow the same correlation: the ratio 
of branching ratios is driven by fB s /fB d \Vts/Vtd\ which 
is determined to a high accuracy in the global fit. On 
the other hand, the NNLO prediction is both lower and 
more accurate than the NLO case, in agreement with the 
results in ref. [25]. This highlights the importance of a 
precise measurement of this observable, e.g., at Bclle-II. 

The study of the time-dependent decay rates of 
B —► I? ± 7r =F , D^tt^ and D ± p T provides a measure 
of rsin(2/3 + 7 ), where r is the ratio of the magni¬ 
tudes of the doubly-Cabibbo-suppressed and Cabibbo- 
favoured amplitudes [38] . Because of the smallness of 
this ratio for the three modes, one has to extract them 
from B° —► D^ + h~ decays assuming SU( 3) flavour 
symmetry (allowing for SU{ 3) breaking at the level of 
1 ± 0.10 ± 0.05). Another additional input needed is the 
ratio of decay constants for excited mesons: /£, /fp = 
1.16 ± 0.02 ± 0.06 [5D[. Combining those observables, we 
obtain a constraint on] sin(2/3 + 7)|, which corresponds to 
a lower limit | sin(2/3 + 7)] > 0.69 at 68 % CL (Fig. |9j). 

IV. NEW PHYSICS IN AF = 2 
A. Additional inputs and parameters 

As discussed in refs. [ 20 l l 22 l l23l I40H46] . neutral-meson 
mixing is a particularly interesting probe of NP. The 
evolution of the B q B q system is described through a 
quantum-mechanical hamiltonian H = M q — iT q /2 as 
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p-value 



FIG. 8. Prediction on the two dileptonic branching ratios 
B(B S — > fj.fi) and B(Bd —► fifi) coming from the global fit 
(without input on dileptonic branching ratios) compared to 
current experimental information (T8] . B(B S — > fifi) is shown 
removing the (1 + y s ) = 1.07 increase due to time integration. 
The NNLO computation in ref. [28] is indicated in colours, 
whereas the NLO computation used in ref. [37] is the region 
delimited by the dashed line. Regions outside the coloured 
areas are excluded at 1 — p > 95.45 %. The oval contours are 
the experimental 1, 2, 3 ... a confidence regions m- 


Observable 

Value and uncertainties 

Ref. 

Asl 

(-47 ± 17) x 1(T 4 

(35] 

a SL 

(1 ± 20) ■ 10" 4 

m 

a SL 

(-48 ± 48) ■ 10~ 4 

ES] 

AF S 

0.081 ±0.008 

□S] 

Bs,B s /Bs,B d 

1.01 ±0.02 ±0.02 

m 

Bs,B s (m b ) 

0.89 ±0.10 ±0.09 

m 


TABLE IV. Experimental and theoretical inputs inputs mod¬ 
ified compared to Ref. [201123] and used in our fits for NP in 
AF = 2. 


the sum of two hermitian “mass” and “decay” matrices. 
Bq — B q (q = d, s) oscillations involve the off-diagonal el¬ 
ements Ml 2 and F\ 2 , respectively. One can fix the three 
physical quantities |M® 2 |, |Ff 2 1 and (j) q = arg(— Ml 2 /Y\ 2 ) 
from the mass difference A M q ~ 2 |Mi 9 | among the eigen¬ 
states, their width difference Ar g ~ 2 | rf 2 | cos (f> q and the 
semileptonic CP asymmetry 

“^ = Im Mt = tSg Sin * = (111 

Resulting from box diagrams with heavy (virtual) parti¬ 
cles, Ml 2 is expected to be especially sensitive to NP [5(7] . 
Therefore the two complex parameters A s and Ad, de¬ 
fined as 

Ml 2 =Mf^ q • A, , A q = |A 9 |e<, q = d, a,(12) 




Y (deg) 


FIG. 9. Combined constraint on 2/3 + 7 using relevant observ¬ 
ables measured in the B —> Dn, D*n (top) and Dp (bottom) 
decays. 


can differ substantially from the SM value A s = Aj_ = 1. 
Importantly, the NP phases <j>^ do not only affect a gj®, 
but also shift the CP phases extracted from the mixing- 
induced CP asymmetries in Bd —> J/ijiK and B s —> 
to 2/3 + (j>j and 2/3 s — <j>£, respectively. There has 
been a lot of interest triggered on this possibility due to 
disagreements with respect to the SM shown first by the 
early measurements from CDF and D0 on the B s mixing 
angle, and further once D0 quoted values of the like-sign 
dimuon asymmetry Asl (measuring a linear combination 
of <Xg L and ds L ). Lat er measurements of the individual 
semileptonic CP asymmetries and mixing angles have 
not been able to explain the D0 measurement, as they 
showed a good agreement with SM expectations. 

In Refs. [20] (233 we h av e determined the preferred 
ranges for A s and A^ in a simultaneous fit to the CKM 
parameters in different generic scenarios in which NP is 
confined to AF = 2 flavour-changing neutral currents. 
We focus here on Scenario I, where we have treated A s 
and Ad independently, corresponding to NP with arbi¬ 
trary flavour structure. In this setting, KK involves 
three other, unrelated, new physics coefficients which will 
not be discussed in the following. We use most of the in- 
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FIG. 10. Complex parameters A^ (up) and A s (down) in 
Scenario I, not including Asl- The coloured areas represent 
regions with 1 — p < 68.3 % for the individual constraints 
(ttexp = a — 4>^/2). The red area shows the region with 
1 — p < 68.3 % for the combined fit, with the two additional 
contours delimiting the regions with 1 — p < 95.45 % and 
l-p< 99.73 %. 


puts involved in the global fit, apart from B(B S —> fifx), 
which is likely to be also affected by New Physics in a 
way that cannot be connected simply to the New Physics 
introduced in A F = 2 processes. In Scenario I, ck is af¬ 
fected by NP independently from the Bd and B s sectors, 
and thus has no impact on the discussion of NP here. 
The remaining parameters can be found in ref. mm- 
One comment is in order concerning the recent re¬ 
assessment of the value of Asl- Members of the D0 


FIG. 11. Complex parameters A^ (up) and A s (down) in 
Scenario I, including Asl. The conventions are the same as 
in Fig. [TOl 

experiment m have considered an additional SM source 
for CP-violation contributing the dimuon charge asym¬ 
metry (coming from the interference of b —¥ ccs decay 
with and without mixing). This modifies the extraction 
of the linear combination of ag L and Og L from the like- 
sign dimuon asymmetry. This correction has been in¬ 
cluded in the latest D0 update, bringing Asl closer to 
its SM value. But the estimate of this correction has been 
challenged later [JHj , as it misses other contributions from 
b —I cus,ucs,uus which could partially compensate this 
new correction. As the theoretical status remains un¬ 
clear, and since Agu has been in the past always very 
difficult to reconcile with the other A F = 2 observables 
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Quantity 

without Asl 

with Asl 

Re(Ad) 

0.94±°;i| 

0 . 88 l° 0 ;?g 

Im(Ad) 

- 0 . 12 tg;“ 

—o.h18;o5 

|Ad| 

o.95±8;i! 

0-89™ 

<t>d [deg] 

—6.9±S;| 

-7-31 il 

Re(A s ) 

i nc+°- 14 
±.UO_o.i3 

1 ni+ 0 - 17 
l-Ul—o.09 

Im(A s ) 

0 QQ+0.04 
U.UO—0.04 

0 02 +0-04 
u.uz _ 0 04 

|A.| 

i nc+°- 14 
-L.UD_o.i3 

1 oi +0 - 17 

i - u± -o.io 

4>s [deg] 


1 q + 2.3 
± ‘°-2.3 

<t>d + 2/3 [deg] (!) 

46-±i2 

38tl° 3 

(j>s - 2/3 S [deg] (!) 

_4Q + 43 
^• — 16 

- 6 ll 43 

Asl [ 10 - 4 ] (!) 


— 7 . 1 I 4.3 

Asl [ 10 ~ 4 ] 

- 

-ioaIj-j 

a SL [ 10 - 4 ] (!) 

— 1T.3±5;| or 1211“ 

—20.71®;® 

«SL [ 10 - 4 ] (!) 

l- 6 ±iJ 

A *°— 1.9 

AIAtps- 1 ] (!) 

0.0028™® 

0.0042l° 0 ;°°f 9 

ATsIps" 1 ] (!) 

o.o9o±g;8S 

0.089l° 0 ° 8 2 2 3 

Ar s [ps _1 ] 

0.08lt°'ogg 

U.Uol _ 0 006 

B-^tu [ 10 “ 4 ] (!) 

0 . 6881 °;”° 

1 0^^+0 065 

i.UOO_o 345 

B —> tv [10 -4 ] 

1.0291°;°°? 

1 0^7+0-062 
l.UO 1 _ 0 . 155 


TABLE V. 68.3% CL intervals for the results of the fits in 
Scenario I, including or not the Asl measurement. The no¬ 
tation (!) means that the fit output represents the indirect 
constraint with the corresponding direct input removed. 


even within our rather generic scenario, we will consider 
two sets of results, with and without the inclusion of the 
D& measurement. 

In addition, we have updated the values of the bag pa¬ 
rameters, following the recent work from the ETMC col¬ 
laboration 35] , working with Nf = 2 dynamical flavours. 
The impact on our study is however small, since their re¬ 
sults showed an excellent compatibility with the previous 
(quenched) study [50] that we used in previous publica¬ 
tions. For the ratio of scalar quenched bag parameters, 
we have assumed that the breakdown between statisti¬ 
cal and systematic errors in ref. [49] followed the same 
pattern as for the SM {B\) bag parameters. All these 
additional inputs are collected in Table m 


B. Constraints on New Physics 

We summarise our results in Tables Ivl and IVII and in 
Figs. EH an dm including or not Asl ■ We find pull values 
for A SL and </>f — 2/3 s of 2.4 a and 2.5 a respectively, 
illustrating the discrepancy between the two constraints 


Quantity(ies) 

Without Asl 
Deviation wrt 
SM Sc. I 

With Asl 
Deviation wrt 
SM Sc. I 

4 >d +2/3 

1.6 a 0.0 0 

1.6 a 0.0 a 

</>? - 2/3, 

0.0 a 1.1 a 

0.0 a 2.6 a 

Asl 

- 

2.7 cr 2.4 a 

«sl 

0.4 a 0.8 a 

0.4 cr 1.1 a 

Osl 

1.0 a 1.0 cr 

1.0 cr 1.0 a 

Ar s 

0.3 a 0.3 a 

0.1 a 0.1 a 

B{B ->• tv) 

1.3 a 0.8 a 

1.3 a 0.2 a 

B(B —> tv), Asl 

- 

2.5 a 2.1a 

4 >f — 2/8 s , Asl 

- 

2.2 a 2.2 a 

B(B ->• tv), ( j>f ~ 2/3 s , A S l 

- 

2.2 a 1.9 a 


TABLE VI. Pull values for selected parameters and observ¬ 
ables in SM and Scenarios I in terms of the number of equiva¬ 
lent standard deviations between the direct measurement and 
the full indirect fit predictions. Two different types of fits, in¬ 
cluding or not Asl are considered. 


in Fig. 11 We do not quote pull values for A md, s in 
Sc. I, as these observables are not constrained once their 
experimental measurement is removed. 

The global constraint on the argument of A s is more 
stringent than what could be assumed by the overlap of 
the constraints from A rrid, A m s and (f>^ — 2/3 s . This 
can be understood as follows: the fit including NP in 
A F = 2 has a discrete ambiguity in the determination of 
p, fj, so that two solutions (symmetrical with respect to 
the origin) are allowed [120, 22]. This translates into two 
possibilities for /3 s b, with opposite signs. The constraint 
from <f>f — 2/3 s also exhibits two preferred solutions for 
arg(A s ). These two solutions cannot be distinguished at 
1 a if only cjrf — 2(3 S is considered, but the degeneracy 
is lifted in favour of the SM-like solutions once the other 
constraints are added, leading to a global constraint cen¬ 
tered around the solution corresponding to the SM-like 
solution for p, fj , with a domain smaller in size than the 
constraint from <j>^ — 2/3 s 

The comparison between the fits with and without 
Asl shows a slight decrease for A,/ when Asl is added, 
whereas |A S | is essentially unchanged. One notices also 
that in the absence of Asl, the predicted value of a SL can 
take two different values (a small negative one or a large 
positive one), corresponding to the two branches allowed 
by + 2/3. The predicted value for </>£ + 2/3 varies sig¬ 
nificantly when Asl is added or not, since it comes from 
the combination of the constraint from a measurements 
with the semileptonic asymmetries. This yields a notice¬ 
able change in the prediction for AT^. Even though the 
predictions for AT S and B —► tv also seem to vary, this 
mainly concerns the best-fit point and is much less the 
case once ler intervals are considered. 

One can also consider the p -value of the SM hypothesis 












12 


following the discussion in ref. [2DI - Without Asl, the p- 
value for the 2D SM hypothesis A a = 1 (A., = 1 ) is 0.9 a 
(0.3 cr), and the 4D SM hypothesis A^ = 1 = A s = 1 has 
a p-value of 0.7 a. With Asl, the p-value for the 2D SM 
hypothesis A^ = 1 (A s = 1 ) is 1.2 a (0.3 cr), and the 4D 
SM hypothesis A^ = 1 = A s = 1 has a p-value of 1.0 a. 

The two complex NP parameters A^ and A s 
(parametrising NP in M^ s ) are not sufficient to absorb 
the discrepancy between the D0 measurement of Asl 
and the rest of the global fit. The situation has how¬ 
ever improved compared to earlier analyses, due to the 
decreased discrepancy of Asl compared to the Standard 
Model. Without Asl, the fit including NP in A F = 2 
is good, but the improvement with respect to the SM is 
limited. In addition, we stress that data still allow size¬ 
able NP contributions in both Bd and B s sectors up to 
30-40% at the 3cr level. 


V. CONCLUSION 

This article collects a selection of SM predictions 
driven by the global fit of the CKM parameters, in view 
of related recent or foreseeable experimental measure¬ 
ments. The results were obtained by combining the in¬ 
puts collected in Table |I| using the statistical frequentist 
framework adopted by the CKMfitter group. The overall 
agreement of the Standard Model global fit is impressive, 
as confirmed by the representation of the various unitar- 


ity triangles and the results given in Table |III| gathering 
the SM predictions using the inputs. We discussed the 
status of some quantities of importance whose status has 
changed recently ( a , 7 , B(B S —>• pp)). We also provided 
predictions for various observables of interest, as well as 
a table of pulls. 

We have also performed a global fit to flavour physics 
data in a scenario with generic New Phiscs in the BdBd 
and B S B S amplitudes, as defined in Refs. HOI nil- The 
discrepancy between A$l and the rest of the neutral- 
meson mixing observables remains even in this extended 
scenario. If we remove Asl> because of the potentially 
large (and unknown) systematics affecting its extraction, 
the fit improves significantly, with a SM-like scenario be¬ 
ing very likely. However, significant contributions from 
NP are still possible at the 3cr level. This is an invita¬ 
tion for more study of these observables with the LHCb 
upgrade and the start of Belle-II, as discussed in the 
prospective exercise of ref. [ 22 ] ■ 
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